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Abstract 

In this paper we study Doob's transform of fractional Brownian 
motion (FBM). It is well known that Doob's transform of standard 
Brownian motion is identical in law with the Ornstein-Uhlenbeck dif- 
fusion defined as the solution of the (stochastic) Langevin equation 
where the driving process is a Brownian motion. It is also known that 
Doob's transform of FBM and the process obtained from the Langevin 
equation with FBM as the driving process are different. However, also 
the first one of these can be described as a solution of a Langevin 
equation but now with some other driving process than FBM. We are 
mainly interested in the properties of this new driving process denoted 
Y' 1 '. We also study the solution of the Langevin equation with 
as the driving process. Moreover, we show that the covariance of 
grows linearly; hence, in this respect yw is more like a standard Brow- 
nian motion than a FBM. In fact, it is proved that a properly scaled 
version of converges weakly to Brownian motion. 

Keywords: fractional Brownian motion, fractional Ornstein-Uhlenbeck 
process, long range dependence, short range dependence, covariance 
kernel, weak convergence 
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1 Introduction 



It is well known that the Ornstein-Uhlenbeck diffusion U = {Ut ; t > 0} can 
be constructed as the unique strong solution of the Langevin SDE 

dU t = -aU t dt+ dB t , (1) 

where a > and B = {Bt : t > 0} is a standard Brownian motion initiated 
from 0. Solution of (pQ) can be expressed as 

U t = e~ at fx + f\ as dB.Y (2) 

where x is the (random) initial value of U. Using partial integration, the 
stochastic integral in ([2]) can be written as 

e au dB u = e as B s - [ ae au B u du. (3) 
J o 

The stationary distribution of U is iV(0, l/2a). Consequently, taking x to 
be a normally distributed random variable with mean and variance l/(2a) 
independent of U gives us a stationary version of the Ornstein-Uhlenbeck 
diffusion. 

Let B^ = {b\ ' : t > 0} be another standard Brownian motion initi- 
ated from and independent of B. Introduce for t E R 



B t 



B t , t > 0, 
B^ , t < 0. 



The process B is sometimes called two-sided Brownian motion through 0. 
It is easily seen that 



J —c 



" - 

oo 



i := / e as dB. 



is a normally distributed random variable with mean and variance l/(2a). 
Since lim^-co B s /s = a.s., it follows via, e.g., ([3]) that ^ is well defined. 
Choosing now x = £ allows us to write the stationary solution of ([T]) in the 
form 



t 

at 



Ut = e~ at e' 



There is also another well known construction of the Ornstein-Uhlenbeck 
diffusion. This is due to Doob [1] and expresses the stationary Ornstein- 
Uhlenbeck diffusion U (with time axis the whole R) as a deterministic time 
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change of a standard Brownian motion: 

U t = e~ at B at , t£R, (4) 

where a > and at '■= e 2at /2a. The covariance of U is easily obtained 
from (HD 

E(U t U s ) = ^-e~ a( ~ t - s \ t>s. (5) 

In this note we study fractional Ornstein-Uhlenbeck processes. These are 
processes constructed as U above but now the Brownian motion is replaced 
with the fractional Brownian motion (FBM). It is known that the process 
obtained as the solution of the Langevin SDE with FBM as the driving 
process does not coincide with the process obtained as Doob's transform 
of FBM. In Cheridito et al. [3j it is proved that the covariance of the 
former one behaves like the covariance of the increment process of FBM. In 
particular, if the Hurst parameter H is bigger than 1/2 the process is long 
range dependent. On the other hand, the covariance of Doob's transform^ of 
FBM decays exponentially and, hence, the process is short range dependent 
for all values of H E (0, 1). Our main contribution in this paper is to extract 
from Doob's transform the driving process, to study its properties and use 
the process in the Langevin SDE to generate new kind of fractional Ornstein- 
Uhlenbeck processes. 

In the next section we discuss the basic properties of FBM important 
for our purposes. To make the paper more readable, we also recall some 
results from [3]. In the main section of the paper the new driving process is 
constructed and the solution of the associated Langevin SDE is introduced. 
The covariance of the driving process and also the covariance of the solution 
have kernel representations in case H > 1/2. It is proved then that the 
driving process and the solution are short range dependent. Moreover, it 
is seen that it is possible to scale the driving process so that it converges 
weakly to a Brownian motion as the scaling parameter tends to infinity. 

2 Preliminaries 

2.1 Fractional Brownian motion 

Let Z = {Zt : t > 0} be a fractional Brownian Motion, FBM, with self- 
similarity (or Hurst) parameter H E (0, 1), that is, Z is a centered Gaussian 

1 In [3] this transform is called Lamperti's transform (see Lamperti [7]). 
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process with the covariance function 

B(Z t Z s )= 1 -(t 2H + s 2H -\t-s\ 2H ). (6) 

Notice that 

E(Z$) = and B(Zf) = 1, 

and, hence, in particular Zq = 0. Using Kolmogorov's continuity criterion 
it can be proved that Z has a continuous version; therefore, we take Z to 
be continuous. In fact, Z is locally Holder continuous of exponent a for all 
a < H. 

Fractional Brownian motion is if-self-similar in the sense 

{Z at : t > 0} = {a H Z t :t>0} for all a > 0, (7) 

where = means that the right hand side and the left hand side are identical 
in law. This follows from ([6]) because the covariance function determines a 
mean zero Gaussian distribution uniquely. Moreover, from ([6]), for t 2 > t\ > 

s 2 > si 

E((Z t2 -Z tl )(Z S2 -Z Sl )) 

= \ ((ta - si) 2 " - (ti - Sl f H - (t 2 - s 2 f H + ( tl - s 2 f H ) . (8) 

Since the function 

s ■-> (t 2 - s) 2H - (h - s) 2H , s<h< t 2 , 

is decreasing for H > 1/2, and increasing for H < 1/2 it follows that the 
increments of Z are 

• positively correlated if H > 1/2, 

• negatively correlated if < 1/2. 

Consider now the increment process of Z defined as 

I z := {Z n+ i - Z n : n = 0, 1, 2, . . . }. 

It is easily seen that Iz is a stationary second order stochastic process and, 
from ©, 

p Iz {n) := E (Z 1 (Z n+1 - Z n )) = H(2H - IJtT 2 ^ + 0{n 2H ^). (9) 
Next we recall the following definition (see Beran pQ p. 6 and 42). 
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Definition 1. Let X = {X n : n = 0, 1, 2, . . . } be a stationary second order 
stochastic process with mean zero and set px( n ) '■= E (XiXi +n ) , where i is 
arbitrary non-negative integer (by stationarity, px( n ) does not depend on i). 
Then X is called 

{i) long range dependent if there exist a E (0, 1) and a constant C > 
such that limn^oo px(n)/(C n~ a ) = 1, 

(ii) short range dependent if lim^oo J2n=o Px{ n ) exists. 

From Definition [I] and formula Q it follows that the increment process Iz 
of the fractional Brownian motion Z is 

• long range dependent if H > 1/2, 

• short range dependent if H < 1/2. 
Notice that, since Zq = 0, we have 

oo 

lim E(Z N Z 1 ) = Y < Pi z (n). (10) 

n=0 

2.2 Fractional Ornstein-Uhlenbeck processes of the first kind 

We replace now the Brownian motion B in (pQ) with the fractional Brownian 
motion Z, and consider the SDE 

du (Z,a) = _ aU {Z,a) dt + ^ (n) 

Analogously with ([2|), the solution can be expressed as 

Ul z ' a) ( z ) = e- at (x + J* e as dz)j (12) 

with some (random) initial value x. The stochastic integral exists pathwise 
as a Riemann-Stiltjes integral (see Cheridito et al. [3]) and it holds 

/ e au dZ u = e as Z s - \ ae au Z u du. (13) 
Jo Jo 

Furthermore, we introduce Z, two-sided fractional Brownian motion through 
0, and consider 

£:= f e as dZ s . (14) 

J — oo 
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Recall that the process given by 

7 {o) Jo, t = 0, 

1 ' \t 2H z 1/t , t>0. 

is again a fractional Brownian motion. Therefore, 

lim Z s /\s\ 2H = a.s., 

and, consequently, £ is well-defined (via (|13p ). Taking in (|12p x = £ we write 
the solution in the form 



(Z,a) _ -at ' ' 



ur^> = e 



! e as dZ s . (15) 

J —oo 



Since the increments of Z are stationary and the stochastic integral is a 
Riemann-Stiltjes integral it follows that the process is stationary. 

The stationary probability distribution, i.e., the distribution of £, is normal 
with mean and variance (see Cheridito et al. [3]) 

F(2H + 1) sinM) _ 2H f +0 ° Ixl 1 - 211 , 
a — dx. 



7r 7 1 + x 2 

In case H = 1/2, the variance equals l/2a, as it should. 

Definition 2. XTie process given in ( fi5|) is called the stationary 

fractional Ornstein-Uhlenbeck process of the first kind. 

Next we recall the asymptotic formula for the covariance of C/( Z < Q ) taken 
from [3J Theorem 2.3, which is then applied to derive the range dependence 
properties of U^ z,a \ 

Proposition 3. Let H G (0, |) U (|, 1] and N = 1,2, Then for fixed 

s£R and t — > oo, 

E(C/f'«)^) 

AT /2n-l \ 

= | E a_2n II ( 2H - k n t2H ~ 2n + o{t 2H ~ 2N - 2 ). (i6) 

n=l \ fc=0 / 

Proposition 4. XYie stationary sequence {XJ^ : n = 1,2,...} (and, 
equivalently, the process U^'^) is long range dependent when H > 1/2, 
and short range dependent when H < 1/2. 
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Proof. Leading term of the sum in (|16p is of the order t 2H 2 . Consequently, 

oo oo 

n=0 n=0 

oo 

* 5> 2 "" 2 , 

n=0 

which, by Definition [H gives the claim. □ 



3 Fractional Ornstein-Uhlenbeck processes of the 
second kind 

3.1 Definition and some basic properties 

In this section we derive from Doob's transform of Z a Gaussian process 
with stationary increments. This process is used as the driving process in 
the Langevin SDE. In this way we construct a new family of Gaussian pro- 
cesses which we call fractional Ornstein-Uhlenbeck processes of the second 
kind. This terminology can be justified by observing that in the standard 
Brownian case, i.e., H = 1/2, these processes coincide with the Ornstein- 
Uhlenbeck diffusions; as also do the fractional Ornstein-Uhlenbeck processes 
of the first kind introduced in Definition [2j 

Doob's transform of Z is the process given by 

X [ D ' a) :=e- at Z at , feR, (17) 

where a > and at '■= a(t,H) := H e at ^ H /a. The covariance of X can be 
computed from ([6]). Indeed, for t > s we have 

E(xl D > a) x( D >°)) (18) 

since is a Gaussian process it follows herefrom that it is stationary. 

In particular, using the self-similarity property of the fractional Brownian 
motion (see ([7])) it is seen that xj: D ' a ^ is for all t normally distributed with 
mean and variance (H/a) 2H . 

Proposition 5. The stationary process {X { t D ' a) : t £ R} is, for all H G 
(0,1), short range dependent. 
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Proof. Formula (|18p yields for a fixed s as t — » oo 

E(X t (D ' Q) Xi D ' Q )) = 0(exp(-a min{l,(l -#)/#}*)), (19) 
and this implies the result. □ 
Consider now the process defined via 

Y t (a) := f\~ as dZ as , (20) 
J o 

where the integral is a (pathwise) Riemann-Stiltjes integral (cf. Section f2 . 2 [) . 
In case H = 1/2, y( Q ) is, for all a, by Levy's theorem a standard Brownian 
motion. Using y( a ) the process can be viewed as the solution of the 

equation 

dx (D,a) = _ aX (D,a) dt + ^(a) (21) 

Proposition 6. For all a > 

{a ff y^ : t > 0} i : t > 0}. (22) 

77ie process y( a ) has stationary increments. 
Proof. Integrating by parts we obtain 

Y t (a) = f e~ as dZ as = e~ at Z at - Z ao + a f e~ as Z as ds (23) 
J o J o 

Using ([?D - the self-similary property of FBM - the claimed identity in law 
([22]) follows from (f23j) . Moreover, the equality 

e {{y t [ a) - y t [ a) ) - ^)) = e ((y t ( 2 i - y&) (ia - 

holds for t% > t\ > S2 > s\ > and h > again by the self similarity of FBM 
and exploiting (|23|) . Consequently, the increments of y( Q ) are stationary. □ 

Inspired by Proposition [6J we consider the Langevin SDE with MS 
the driving process: 

dui Dn) = - 7 U t {Dr,) dt + dY t {1 \ 7 >0. (24) 



S 



The solution can be expressed (cf. (fT5|) ) 

J — oo 

= e~T* T e^- 1 )'^, 7 >0, (25) 

./ — oo 

where K^ 1 ) stands for the two sided Y^ process. To show that the stochastic 
integral term makes sense also for 7 £ (0, 1] recall first that for all (5 < H 

KmZ 8 /\sf = a.s. (26) 

because Z is Holder continuous of order (3 < H. Next for T < using partial 
integration 

f e^" 1 )" dZ au = e^-V'Za, - e^ T Z aT - (7 - 1) f e^ u Z au du, 
Jt Jt 

and by (|26p the right hand side has a well defined limit as T —* — 00. 

Since the increment process of is stationary it follows that U^- D '^ is 
stationary and, therefore, we have well justified the following 

Definition 7. The process U^ D '^ defined in Ii25\) or, equivalently, via the 
SDE \21$ is called the fractional Ornstein-Uhlenbeck process of the second 
kind. 

We conclude this section by characterizing the Holder continuity of y( a ) 
and U^ D,1 \ The result holds for more general stochastic integrals with re- 
spect to Z (see Zahle 0), but the following simple proof in our special case 
is perhaps worthwhile to present here. 

Proposition 8. The sample paths ofY^ and U^ D '^ are (locally) Holder 
continuous of order j3 < H. 

Proof. Prom (|2ip we have 

Y (a) = x (D,a) _ x {D,a) + j* aX ( D , a) ^ (2?) 

J 

Consequently, t 1— > is continuous and the Holder continuity properties 
of y( Q ) and X( D <°) are the same. Hence, let T > be given and consider 
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for s, t < T and (3 > 



x (D,a) _ x (D,a) 



\e- at Z, 



at 



t-s 



-at 



Za s \ + Z as e" 



-at g— as 



< 



Za s \ \<H ~ a. 



t - s 







\a t - a s 



\7 —7 

< Kt-—^— as 



t-sf 
+ C T , 



+ |t- s | 1 "' 3 |Z as 



| e — at _ g — as | 

It - sl 



I at - dsf 

where and Ct are (random) constants which do not depend on s and 
t. The claim follows now from the fact that the paths of FBM are (locally) 
Holder continuous of order (3 < H. Similarly, for the process U^ D ^ (starting 
from 0) we have 

rt 



Y t (1) -7e" 7 * f e^Y^ds, 
Jo 



and it follows that also U^ '^ is Holder continuous of order [3 < H. □ 

3.2 Kernel representations of covariances and short range 
dependence 

We make now the following assumption valid throughout the rest of the 
paper 

1/2 < H < 1. 

In this case, as is easily checked, the covariance of the fractional Brownian 
motion has for t<i > t\ and S2 > s\ the kernel representation 

f*2 rs2 



E ((Z t2 - Z tl )(Z S2 -Z S1 ))= I" f 2 H(2H - l)\u - v\ 

Jt\ J Sl 



2H-2 



du dv. 



In the next proposition we derive an analogous representation for the process 
The result is formulated for all values on a > 0. 

Proposition 9. The covariance ofY^ with 1/2 < H < 1 has the kernel 
representation 



E ((yW_yW) (yw-yw)) 

rt2 I-S2 



-a(l-H){u-v)/H 



dudv, (28) 
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where > t\, S2 > s\, and 

C(a,H):=H(2H -!)[-) 

The kernel 

e -a(l-H)(u-v)/H 
r a , H (u, v) := C(a, H) - _ e _ a(u _„ )/H|2(1 _ g) (29) 

is symmetric, i.e., r a H(u,v) = r a) u{v,u) for all u,v G R. 

Proof. Recall the formula (see Gripenberg and Norros [5] Proposition 2.2) 

= H(2H-1) f [ f(s)g(t)\s-t\ 2H - 2 dtds, (30) 
Jn. Jn 

where 1/2 < if < 1 and /, g e L 2 (R) n L X (R). Since 

n W == / e" as ^ as = e" at Z at - Z ao + a I e~ as Z as ds 
Jo Jo 

simple manipulations yield 

Y W=H H s- H dZ s . 

J CLQ 

Consequently, the claim follows by a straightforward application of (|3U|) . □ 

Remark 10. Notice that the kernel r Q) # is in L 2 ([0,T] x [0, T\) if and 
only if H > 3/4. Consequently, forY^ we have similar absolute continuity 
properties as for fractional Brownian motion (see Cheridito ]2$)- Namely, 
the measure induced by the process {B t + y} 1] :t>0}, where y (1) and the 
Brownian motion B are assumed to be independent, is absolutely continuous 
with respect to the Wiener measure. 

For the next result, recall from Proposition [6] that the increments of 
are stationary. 

Corollary 11. The increments ofY^ are positively correlated. The incre- 
ment process ly := {i^+i — Yn\ n = 0, 1, . . .} is stationary and short range 
dependent. 
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Proof. From (|28p it follows immediately that the increments are positively 
correlated. Of course, we may also deduce from (|28|) the stationarity of the 
increments of 7", To show that ly is short range dependent consider 



rn+l rl 

E (1^(1^-1 - K)) = du dv r a>H {u, v) 

Jn J 



-a(l-H)n/H 

du I dve -^H)(u-vy Hll _ e -n/H e -a(u-v)/Hm-l)_ 



The integral term has a positive finite limit asn-> oo. Indeed, Lebesgue's 
dominated convergence theorem yields 

lim f l du f 1 dve -^-H)(u-v)/H ][l _ e -n/H e - a {u-v)/H ] 2{H-l) 
n-*oo J J 

"'du I' dve-< l - H ^-^ H . 



Consequently, 

PYia) (n) := E (yW(i£l - Y^)) = O (e^^W") . (31) 
and, hence, 



lim E (4 a) y/ Q) ) = PyW in) < +oo (32) 

n=0 

completing the proof. □ 

Next we study the asymptotic behaviour of the variance and covariance 
of Y^ a \ For this, it is practical to rewrite the symmetric kernel r aj n m d2H]) 
as 

r a ,H(t,s) = k a ^ H (t - s) 

with 

k a>H {x) := C(a,H)e- a ^- H ^ H \1 - e~ ax/H \ 2H - 2 . (33) 



Proposition 12. The following formulas hold: 

E((y/ Q) -yi Q )) 2 ) =2 J* \t-s-x)k a , H (x)dx, (34) 
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E (Y t (a) Yj a )) = _/ (* - *) *a,ff (*) dx (35) 

/■s rt—S 

+ / (s — x) k a ^H (x) dx — / {t — s — x)k a ^H{x) dx. 
Jo Jo 

Moreover, 

E f(y t (a) ) 2 ) = 0(t) as i -» oo, (36) 



and 

lim E(y/ a) y s (a) ) = s I k a H (x)dx + f {s-x)k a H {x)dx. (37) 
' Jo Jo 

Proof. We apply fl28|) to obtain (|53]) : 

t ft 



du J dvr a , H {u,v) 
du / dv r a ^u{u,v) 

J s 

ft ry-s 
= 2 dy dx k a ,H(x) 

J s Jo 

rt—s rt 

= 2 / dx I dy k a>H (x) 

Jo Jx+s 
rt—s 

= 2 / [t — s — x) k a> H{x) dx. 
Jo 

Putting here s = and using 

/■OO /'OO 

/ k a ^n{x)dx <oo and / x k a> H(%) dx < oo 
io Jo 



yield (|36p . Furthermore, straightforward computations produces formula 
(|35p from (|34p . It remains to to prove (|37p . Consider for t > 2s 

E((y/ a) -yi«))yi Q )) = /"'du f dvr a , H {u,v) 



s JO 
t PS 

du / aY> k a: i{(u — v) 
Jo 

rt—s 

x k a h (x) dx + s / k a ^H(x)dx 







+ / (t — x) k a ,H{x) dx. 

t-s 
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Consequently, 

limE ((yW-yW)yW) 



rs poo 

/ x k aj H(x) dx + S / k a ^ff(x) 

JO Js 



DO 

dx 



OO PS 



t— >oo 



= s k a> H(x)dx— / (s — x) k a} H(x) dx, 
Jo Jo 

from which (|37p easily follows. □ 

Remark 13. The short range dependence property ofY^ also follows from 
( ~37\ ) since (recall that Y Q (a) = 0) 



f> yCa) (n) = Jim E (Y^Y} a) ) < +00. 



n=0 



Proposition 14. The covariance of U^ D '^ has the kernel representation 

, t r s p(7-i+^)(«+«) 



00 J —00 



Proof. As in the proof of Proposition [9j we use also here formula (|30p . 
However, now we need an extended version due to Pipiras and Taqqu [8] 
stating that (i30j) holds true for functions / and g satisfying 



\f(s)\\g(t)\\s -t\ 2H ~ 2 dtds < 00. (38) 
Consider 

/t pat 
^- l >dZ as = H-^ H s^ H dZ s . (39) 
-00 Jo 

To check that condition ([38]) is valid for f(s) = g(s) = s( 7_1 ^l( 0i ai)(s) it is 
enough to show that 

11 uvp-^ H \u-v\ 2H - 2 dudv 



JO 



2 I' duu^ H I dvv^ H (u-v) 2H - 2 <cx 



JO 
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The inner integral can be expressed in terms of the Beta-function 

[ U ( u _ v fH-2 dv = u (7+i)H-i Be ta(l + (7 - 1)H, 2H-1). 

Jo 



Consequently, 
-1 rl 



! ! (uv) { ^ H \u-v\ 2H - 2 dudv = — Beta(l + (7 - l)il, 2# - 1) 
Jo Jo lH 



and condition (j38[) holds. To verify the claimed kernel representation is now 
a straightforward computation using formula (|30p . □ 



Recall from Corollarvl 1 1 1 that the increment process of Y"W is short range 
dependent, and that if is used as the driving process in the Langevin 
equation the solution is the process U^ D '^ . In the next proposition we show 
that also U^'^ is short range dependent. Formula (140p can be compared 
with the corresponding formula (|19p for X^ D ' a h In fact, (|19p with a = 1 is 
(f40j) with 7 = 1, as it should. 



Proposition 15. The rate of decay of the covariance ofU^ D '^ is exponen- 
tial. More precisely, 

E (V t (D ' 7) [/ s (D ' 7) ) = O (exp (- min{ 7 , (1 - H)/H} t)) , as t — ► 00. (40) 

In particular, the stationary process range dependent. 

Proof. Without loss of generality, we may take s = and, hence, consider 



00 ^ —00 



l e u/.ff _ e i)/H|2(l-H) 



= Ai(t) + A 2 (t), 
where, for some fixed T > 0, 

,T ,0 P (7-l+ W )(«+«) 



L (t) : =H 2H - 2 e-^ f du f dv 
J —00 J— 00 



| e «/H _ e t)/H|2(l-H) 



and 



A 2 (t) := H 2H - 2 e~~ ft f du f dv 
Jt J -co 

Clearly, 



t f0 e (7-l+^)("+«) 



Ai(t) = 0(exp(-7i)) as t -> +00. 
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For the integral term in A2(t) we have 



t r o e ( 7 _i + ^)(«+„) f . t r o P (i+i-jj)u p ('r-i+jj)v 



T d - y_ ^ | e ^-e^pd-H) = y T d - (1 _ e( „_ M)/ H )2 (i-H) • 

For G (T,t) x (-oo,0) 

1 < M - e ( ^ u)// M < f 1 - e- T/H J 

and, consequently, formula ((ID]) holds. □ 

3.3 Weak convergence of to Brownian motion 

In Proposition [12] it is proved that the growth of the variance of is 
asymtotically linear as t — > +oo (see ([36]) ). This suggests that yM, when 
properly scaled, behaves asymtotically as a standard Brownian motion. We 
give the precise statement in the next proposition formulated for arbitrary 
a > 0. 

Proposition 16. For a > define 

\/a 



and let B = {B t : t > 0} denote standard Brownian motion started from 0. 
Then as a — > +oo 

{Z, (a ' Q) : i > 0} we 4 kly : t > 0}, 

where w =£ ly stands for weak convergence in the space of continuous functions 
and a = o~(a,H) is a non-random quantity depending only on a and H (see 

Proof. We show first that the finite dimensional distributions of Z^ a,ct ' con- 
verge to the finite dimensional distributions of o~B. Since is a Gaussian 
process with mean zero it is enough to verify the convergence of the covari- 
ance function. From (]35]) in Proposition [12] we have for t > s 



a 
1 

a 



at 







(at — x) k a> H 0*0 dx + / (as — x) k a ^H(x) dx 
Jo 

ra(t—s) 

(at — as — x)k a ^H(x) dx ) 
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with k a> H defined in ([33]) . Letting here a — > +00 yields, after some simple 
computations, 

poo 

lim E Zrz{ a) ) = 2s k nH (x)dx 



j t £jl'\ = Z S / Kq, H 

a— >oo \ / 

= k(ch,H)s, 

where 

As(a, fl") := 2 C(a, £T) — Beta(l - H,2H - 1) 

a 

and Beta(l — H, 2H — 1) is the Beta function. Since E(B t B s ) = s for t > s 
we have proved the convergence of finite dimensional distributions of Z^ a '°^ 
to the finite dimensional distributions of aB with 



a = a(a, H) = \J n(a, H). (41) 

To prove tightness, it is enough to verify (see, e.g., Lamperti [6]) that there 
exists a constant C (might depend on a and H) such that for all a > and 

t > s 



A := E ((Z t (a ' a) - Z^Y^j <C 



(t-s). 



We have by formula (|34|) in Proposition 

2 



A = iE((yW-ywy 



1 Mt-s) 
2— / (a(t — s) — x) k a ^H {x) dx 



a 

<C(t-s) 



with, e.g., C = 2 L°° k aj ff(x) dx. This completes the proof. □ 
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